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1 Introduction 

The refinement equation is a functional equation of the form 

N 

f{x)=Y,c j f{Xx-d j ), (1) 

3=0 

where A > 1 and all the Cj,dj are real numbers. For the refinement equation 
((T|), the value A is called a dilation, whereas the numbers {dj} are referred to as 
translations. Throughout this paper, we suppose that do < d\ < • ■ • < djv and 
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define the Fourier transform of f(x) by the formulae 

/oo 
f(x)e~ 2 ^ x dx. 
-oo 

Taking the Fourier transform of both sides of (HJ) we obtain 

f(w) = H(\~ 1 w)f(\~ 1 w), (2) 

where H(w) = A -1 J2f=o Cje~ 2lTldjW is called the mask polynomial of the re- 
finement equation. Setting w = into l[2]). we obtain X] ^Lo c i = ^ provided 
that /(0) 7^ 0. For simplicity, we shall call the function f(x) satisfying (JTJ) with 
/(0) = 1 a X-refinable function with translations {dj | < j < iV}. It plays a 
fundamental role in the construction of compactly supported wavelets and in 
the study of subdivision schemes in CAGD ([2[7])- 

The existence and regularity of the refinable function are of some interest. 
These questions were studied on several occasions. In [6 , Daubechies and La- 
garias showed that up to a scalar multiple the refinement equation ([1]) has a 
unique distribution solution / satisfying supp/ C [do(\ — l) _1 ,djv(A — 1) ]. 
Moreover, they also showed the nonexistence of a C°°-refinable function with 
compact support in one dimension when the number A and all the dj are in- 
tegers. In ;2], Cavaretta et al. extended this result to higher dimensions by a 
matrix method. When A is non-integer, 'the regularity question becomes more 
complicated and perhaps more interesting from the viewpoint of pure analysis' 
[2]. Moreover, the refinable functions with non-integer dilations play an impor- 
tant role in the construction of wavelets with non- integer dilations [T]. Hence, it 
has attracted a considerable attention. For example, the regularity of Bernoulli 
convolutions, which are solutions to 

f(x) = ±f(\x) + ±f(\x - 1), 

was already studied in p~21 [TU [22J, [25] . In general, one characterizes the regu- 
larity of f(x) by considering the decay of f(w). In [J], Dai et al. considered 
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the uniform decay of the Fourier transform of the refinable functions with non- 
integer A. In particular, they give an elegant answer to the following question: 
for any given dilation factor A > 1 and a positive integer k, can one construct 
a X-refinable function f e C k ? In this paper, we reverse the question: can 
one find a dilation factor A > 1 such that there exists a compactly supported 
X-refinable function f € C°° ? 

This question is interesting for two reasons. First, when A and all the dj 
are integers, the regularity of refinable functions is related to the spectrum of a 
certain matrix that is constructed by the refinement coefficients Cj. In this case, 
a negative answer to this question was given in [5J. However, the matrix can 
only be constructed in the case when A and all the dj are integers. Therefore, 
it is not clear whether there exists a non-integer refinable function which is 
compactly supported and belongs to C°°. Second, as pointed out in Section 
3, the regularity of refinable functions is closely related to the distribution of 
powers of a fixed number modulo 1, which is a classical problem in number 
theory. Hence, overlooking the regularity of refinable functions with non-integer 
dilations one will miss out on the beautiful connection between analysis and 
number theory. 

In this paper, combining the tools of number theory with some results of 
harmonic analysis, we extend the result of Daubechies and Lagarias to the 
general case, thus giving a negative answer to the above question when all the 
dj are rational numbers and A > 1 is a real number. 

Theorem 1.1. Let dj G Q and A > 1. Then the refinement equation 

N 

f(x) = ^2cjf(Xx-dj), 

3=0 

has only the trivial compactly supported C°° solution, i.e., f = 0. 

The proof of Theorem 11.11 is based on a purely number-theoretic statement 
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concerning the distribution of powers of a fixed number modulo 1 (see Theo- 
rem 13.11 below) . 

Remark 1. It should be noted that any translate of a refinable junction is 
refinable. If f(x) satisfies (QJ) theng(x) = f(x—b/(X—l)) satisfies the refinement 
equation 

N 

g( x ) = c i9( Xx - d 3 + &)< 

3=0 

which has the same dilation, but a different translation set {dj — b | < j < N}. 

In Section 2, using Theorem ll.il we extend the results of [3] concerning the 
refinable splines with non-integer dilations. We also construct a counterexam- 
ple to the conjecture about the refinable spline with non-integer translations. 
The box spline B(x\M) associated with the s x n matrix M is considered as 
multivariate generalization of the univariate B spline and have become 'one of 
the most dramatic success of multivariate splines' |10j . It is well-known that 
the box spline is a A-refinable function for any positive integer A > 1. Hence, 
box splines play an important role in the subdivision algorithm and wavelets. 
Moreover, one is also interested in the question when a refinable distribution is 
essentially a box spline [15l [26] . However, very little is known about box splines 
satisfying the refinement equations with non-integer dilations. In Section 2, we 
give a characterization of refinable box splines with non-integer dilations. 

The proof of Theorem 1 1 . 1 1 will be given in Section 3. Other proofs concerning 
refinable splines will be given in Sections 4 and 5. 

2 Refinable Spline Functions 

In this section, we study one of the most interesting refinable functions, so-called 
spline function. Firstly, let us recall the definition of spline functions, B splines 
and box splines. 
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Definition 2.1. The spline function f(x) is a piecewise polynomial function. 
More precisely, there are some points — oo = Xq < X\ < ■ ■ ■ < % < ^Af+i = 
+00 and polynomials Pj{x) such that f(x) = Pj(x) for ) for each 

j = l,...,M+l. 

We call the points Xj, j = 1, . . . , M, the knots of f(x) and maxj{deg(P,)} 
the degree of f(x). Splines are widely used in the approximation theory and in 
computer aided geometric design. The spline functions satisfying a refinement 
equation are the most useful ones. A special class of refinable splines are so- 
called cardinal B splines which are defined by induction as 



1 if x € [0,1), 
otherwise, 



B Q (x) = 
and for k > 1 

Bk = Bk-i * B , 

where * denotes the operation of convolution. By the definition of cardinal B 



splines, the Fourier transform of B splines is given by the formulae 
B k {w) 



1 — exp(— 2niw)^ k 



2niw 

A simple observation is that Bq is m-rcfinable for any integer m > 1 and 
satisfies 

m— 1 

B o{x) = B (mx-j). 
3=0 

Since Bk = B * B a * ■ ■ ■ * B , the B spline B]~ is also m-refinablc for any integer 
m > 1. 

The box spline B(x\M) is defined using Fourier transform by the formulae 

S^.gizf^pa), (e c; (3, 

where M = (mi, . . . , m n ) is an s x n real matrix of (full) rank s. In particular, 
for M = (1, . . . , 1) € Z fe+1 , the box spline B{x\M) is reduced to the univariate 
B spline Bk- 
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Next, we turn to the general refinable spline function. In [5T], Lawton et 
al. gave the characterization of compactly supported refinable univariate splines 
f(x) with the additional assumption that the dilation factor is an integer and 
all the translations are integers. They proved the following theorem: 

Theorem 2.2. f [21]) Suppose that f(x) is a compactly supported spline function 
of degree d. Then f(x) satisfies the refinement equation 

N 

f{x) — Cjf(mx — dj), m > 1, dj G Z, 

3=0 

if and only if f(x) = ^2n =0 p n Bd(x — n — <io/( m — 1)) f or some K > and {p n } 
such that the polynomial Q(z) = (z — l) d+1 ^2^ =0 p n z n satisfies Q(z)\Q(z m ). 

Here and below, Q(z)\Q(z m ) means Q(z) divides Q(z m ), namely, that the 
quotient Q(z m )/Q(z) is a polynomial. (In case the notation a\b is used for 
integers a and b, it means that the quotient b/a is an integer.) 

Remark 2. Theorem \2.2\ was extended to higher dimensions by Sun 1261. More- 
over, in fHy, the authors also gave a characterization of compactly supported 
univariate refinable splines whose shifts form a Riesz sequence. This was gen- 
eralized to higher dimensions in \ 16}/ . In f^, the authors gave complete charac- 
terization of the structure of refinable splines. In \ 15}/ . Goodman gave a review 
on refinable splines (including refinable vector splines). 

However, as pointed out in [3], refinable splines do not have to have integer 
dilations or integer translations. In this case, Dai et al. proved the following 
theorem: 

Theorem 2.3. ( 3 ) Suppose that f(x) is a compactly supported spline satisfying 
the refinement equation 

N N 
3=0 3=0 
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such that A £ R and dj € Z. XTien 

There exists an integer I > such that X 1 € Z. 
f-Bj Lei A; &e the smallest positive integer such that X k E Z. TTien the com- 
pactly supported distribution solution <j)(x) of the refinement equation 

n 

= ^ Cj A fe "V(A fc ^ - dj) (5) 
i=o 

is a spline. 

(C) There exists a constant a such that the spline f(x) is expressible as 

f(x) = a<j>(x) * 4>(\~ 1 x) (piX-^-^x), (6) 

where <j>(x) is the spline given by {5|). 

Conversely, if the refinement equation ^ satisfies (A) and (B) then the 
compactly supported distribution solution f is a spline given in (0|). 

By Theorem 11.11 we can prove Theorem 12.31 under weaker conditions. A 
compactly supported function on R is piecewise smooth if there exist an integer 
M and some real numbers a\ < a 2 < ■ ■ ■ < ajvf such that / £ C°°(oj, aj+i) for 
.7 = 1,2,..., M — 1 and supp/ C [<ii , a M ] ■ 

Theorem 2.4. Let f be a piecewise smooth function with compact support sat- 
isfying the refinement equation 

N N 
3=0 j=0 

where X > 1 and dj € Z. Then f{x) is a spline function. Hence, (A), (B) and 
(C) m TheoremlEMhold. 

To deal with non-integer translations, in [3J, the authors raised the following 
conjecture. Suppose that f{x) is a X-refinable spline that is X-indecomposable. 
Then the translation set for f must be contained in a lattice, i.e., a set of 
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the form aZ + b for some a / 0. Here, we call a A-refinable spline f(x) X- 
indecompo sable if it cannot be written as the convolution of two A-refinable 
splines. As stated in [3J, if this conjecture is true, then one can classify all 
refinable splines by Theorem 12.31 However, we construct a counterexample to 
the conjecture. 

Consider the spline function f(x) = B(x\(l, y/5/2)). Using the Fourier 
transform 

-i p — 2-niw i p —2iri\/b/2w 

f! w ) = — - , 

2niw 2iri^5/2w 

we have H(w) = f(Vl0w)/f(w) = J2j=o e~ 2mdjW , where the numbers 

d ,...,d 9 are equal to 0, 1,2, 3, 4, 5/VlO, 1 + 5/VTO, 2 + 5/VlO, 3 + 5/-y/T0,4 + 

5/VlO, respectively. Hence, the spline f{x) is -v/lO-rcfinable and satisfies the 

refinement equation 

i 9 

with the translation set {0, 1, 2, 3, 4, 5/\/l0, 1 + 5/VlO, 2 + 5/^10, 3 + 5/VlO,4 + 
5/VlO}, which is not contained in a lattice. Moreover, we have the following 
proposition: 

Proposition 2.5. The univariate box spline B(x\(l, y/h/2)) is \f\Q-indecom- 
posable. 

The proof of this proposition is non-trivial, ft is postponed to Section 4. 

Motivated by this counterexample, we shall study the refinable Box spline 
function with non-integer dilation and translations. It is not only helpful in the 
study of general refinable splines with non-integer translations, but also useful 
in understanding of the box splines associated with non-integer matrixes. 

It is well-known that the box spline B(x\M), with an additional assumption 
that M is an s x n integer matrix, satisfies some higher dimension refinement 
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equation with an integer m > 1, i.e., 

B(x\M) = c ( - n) B(mx - j\M), (8) 

where c^™' = m s_ "#{a <G Z n \Ma = mj} (see [9j). For this reason, box splines 
are widely used in computer aided geometry design and wavelets. 

The box splines associated with the non-integer matrix are also of interest. In 
[T51|57], the authors discuss linearly independent integral lattice translates of the 
box splines associated with a non-integer matrix. Here, we shall characterize the 
box splines associated to a non-integer matrix satisfying the refinement equation 
with some non-integer dilation and non-integer translations. 

Theorem 2.6. Suppose that the univariate box spline B(x\A) satisfies the re- 
finement equation 

N 

B(x\A) = J2c j B(Xx - dj\A), 

3=0 

where A > 1, A ^ Z and dj £ R. Then 

(A) There exists a positive integer I such that X 1 £ Z. 

(B) For each element mo £ A, there exists an element m £ A and an integer 
p such that m — pmo/X. 

( C) The vector A contains a sub-vector of the form 

(mo, m pxX~ x , . . . , ra p;_i A~' +1 ), 

where I is a positive integer such that X 1 £ Z, andpj,j = 1, ... ,1—1, are integers 
satisfying pAp J+1 , j = 1, 2, . . . , I - 2, and pi-i\X l . 

(D) If all the dj £ Z, then the vector A can be represented as the union of 
the vectors of the form (mo, Xmo, ■ ■ ■ , A fc_1 mo), where k is the smallest positive 
integer for which X k £ Z and mo £ Z \ 0. 

Moreover, using the technique of dimensional reduction, we can extend the 
above results to higher dimensions. 
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Theorem 2.7. Suppose that the s-variable box spline B(x\M) satisfies the re- 
finement equation 

N 

B(x\M) = CjB(\x - dj\M), 
3=0 

where A > 1 , A ^ Z and dj £ R s . Then 

(A) There exists a positive integer I such that X 1 € Z. 

(B) For each column m in M , there exists a column m in M and an integer 
p such that m = pmo/X. 

( C) The matrix M contains a sub-matrix of the form 

(m ,TOoPiA~\ . . . ,TOop;_iA- (+1 ) e R sxi , 

where I is a positive integer such that X 1 € Z, and pj, j — 1, ... ,1—1, are integers 
satisfying Pj\pj+i, j = 1, 2, . . . , I - 2, and pi-i\X l . 

(D) If all the dj G Z s , then the matrix M can be represented as a disjoint 
union of matrices of the form (m , Atoo, • • ■ , A fe_1 TO ) € R' sx ' c ; where k is the 
smallest positive integer such that X k € Z and mo & Z s \ 0. 

We conclude this section by stating a conjecture which classifies all refinablc 
splines in terms of box splines. A function P(w) is a real quasi-trigonometric 
polynomial if it is of the form of X^jLo c i e ~ 2ntdjW , where Cj,dj € R. The real 
quasi-trigonometric polynomial P(w) is X-closed if P(Xw)/P(w) is also a real 
quasi-trigonometric polynomial. 

Conjecture 2.8. The spline f(x) satisfies the refinement equation 

N N 
f( X ) = C 3-f( Xx - d o)' C 3 = A ' 

if and only if f(w), i.e., the Fourier transform of f{x), can be expressed in 
terms of box splines as follows: 

f(w) = e 2mw P(w)B{w\A), 
10 



where a is some constant, B(w\A) is the Fourier transform of a X-refinable 
box spline B(x\A), and P(w) is a X-closed quasi-trigonometric polynomial with 
P(0) = 1. 

Remark 3. Throughout the paper, we only consider the case where A > 0. 
However, as pointed out in if fix) is a X-refinable spline then it is also a 
(—X)-refinable spline. Hence, the case, where X < 0, can be also studied using 
the results of this paper. 

3 Proof of Theorem 11.11 

Before giving the proof of the theorem, we need a theorem of purely number- 
theoretic nature which plays an important role in the proofs below. For any 
x € R, let ||x||z denote the distance from x to the nearest integer in Z. 

Theorem 3.1. Suppose that X > 1 and n, . . . ,r m are real numbers. Then there 
exists a positive number £ and a positive number c = c(A, r±, . . . , r m ) such that 

||£A'* - T-iHz > c 

for every i = 1, . . . , m and every j = 0, 1, 2, . . .. 

Proof. For each r > 0, put 

S T := U kez U™ ! (k + r t -T,k + n + t). 

Then \\x — r,-||z > t for every i = 1,. ..,mif and only if x € R\S' r . This implies 
that ||£A™ — Ti||z > c (for every i = 1, . . . , m and every integer n > 0) if and 
only if £ £ S c A" n for n > 0. 

We will construct a sequence of closed nested intervals Iq D I \ D Ja D . . . 
(where I C (0, 1)) of length |7 M | = V2A^/A 9M , where M = 0,1,2,..., such 
that £A" ^ S c for every £ e Jm and every n = — 1,0, . . . ,<?M — 1. (Here, for 
convenience, we start with n — — 1, so the final result holds for every i = 1, . . . , m 
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and every n = — 1, 0, 1, ) Also, here g = g(X, m) is the least positive integer 

satisfying 



Then, for the common point £ £ ^m=o^ m ' ^ nc inequalities ||£A™ — r i 1 1 z > c will 
be satisfied for each i = f , . . . , m and each integer n > — 1. 

The proof is by induction on M. We begin with M = 0. Evidently, S c (~\ (0, 1) 
is a union of at most m + 1 intervals of total length 2cm. Thus 5 C A n (0, 1) is a 
union of < m + 1 intervals of total length 2Xcm. The set (0, 1) \ 5 C A is thus a 
union of at most m + 2 intervals whose lengths sum to a number > 1 — 2Xcm. 
It contains a closed interval of length \f2~Xc if 1 — 2Acto > (m + 2)\/2Ac. So the 
required closed interval 7 C (0,1) of length V2Ac exists if 



This inequality clearly holds if c is less than a certain constant depending on A 
and to only. Let us start with this 1q. 

For the induction step M i— > M + 1, we assume that there exist closed 
intervals 7m C Im-i C • • • C 7o with required lengths such that 1m C (0, 1) \ 
S'cA - ™ for each 77 = —1,0,..., 3M — 1. We need to show that Im contains a 
subinterval I M +i of length V2Xc/X3 ( M+1 ) such that J M +i C (0, 1) \ ^A"™ for 
each 77 = gM,gM + 1, . . . ,#(M + 1) - 1. 

Fix 77 e {ffM.^M+l, . . . ,5(M+1)-1}. At most m(l + A™|/ M |) points of the 
form (k + ri)X~ n , where k £ Z and i = 1, . . . , m, lie in Im- So the intersection of 
S'cA - ™ and Im consists of at most 777(1 + X u \Im |) open intervals of length 2c/A" 
(or less) each plus at most two intervals of length c/A" (or less) each at both 
ends of Im- As n runs through gM, . . . , g(M + 1) — 1 the total length of such 
intervals is at most 

S (M+1)-1 ff (M+l)-l 



Im= (m(l+X n \I M \)2cX- n +2cX- n ) = ^ (2cm\I M \+2c(m+l)X- n ). 

n=gM n=gM 
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X a > 2(l+ 5 777). 




< 1. 



U^ng ESm )_1 A- < £~ gM A"* = A^/CA-l) = A|/ M |/(^(A-1)), 
we find that 

f M < 2cTO 5 |/ M |+2c(TO+l)A|/Af|/(V2A^(A-l)) = |/ M |(2cm 5 +V2A^(m+l)/(A-l)). 

The remaining part in Im is of length at least 
\lM\-\lM\(2cmg+V2X^{m+l)/(X-l)) = \I M \{l-2cmg-V2\e{m+l)/{\-l)). 
It consists of at most 

S (M+1)-1 

1+ ^ m(l+A"|/Af|) < l+gm+|/M|A 9 < M+1 V(A-l) - l+.gm+\/2A^A 9 /(A-l) 

n=gM 

closed (possibly degenerate [it, it]) intervals. In order to show that one of these 
closed intervals is of length at least v2Ac/A s( - M+1 ' (so that we can take it as 
Im+i) we need to check that 

|J M |(1 - 2cmg - V2\c{m + 1)/(A - 1)) > \hi+i\(l + gm + V2\c\ 9 /{\ - 1)) 

for c small enough. Indeed, using |/m+i|/|^m| = A -9 , we can rewrite this 
inequality as 

2cto3 + (1 + gm)\- g + V2Xc(m + 2)/(A - 1) < 1. 

Since A 9 > 2(1 + gm), the required inequality would follow from 

2cto 5 + V2Ac"(m + 2)/(A - 1) < 1/2. 

Clearly, g depends on A and m only. So this inequality holds for some positive 
c depending on A and m only. This completes the proof of the theorem. □ 

Remark 4. The m = 1 case of this theorem with some explicit constant c was 
recently obtained by the first named author in Ulf . The existence of such a 
positive number c := c(A) for m = 1 was conjectured by Erdds 17 3f and then 
proved independently by de Mathan JEj and Pollington [23]. In fact, this result 
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was proved already by Khintchine in 1926 (see Hilfssatz III in [191), but then 
forgotten. 

Remark 5. By estimating g from above and by some standard calculations, one 
can see that the inequalities 2Xcm+ (m+2)v2Ac < 1 (corresponding to the M = 
case) and 2cmg + \J 2 Ac(m + 2)/( A — 1) < 1/2 (corresponding to the induction 
step M h-» M + 1 case) both hold if, for instance, c := (A — l) 2 /(20(rn + 2) 2 A 3 ). 
This gives an explicit expression for c in Theorem \3.1\ for each A > 1 . The 
main part for the "difficult" case when A is close to 1 is the factor (A — l) 2 . It 
is essentially the same factor as that in the proof of a similar result obtained by 
the first named author in the m = 1 case \llf - 

We now can give the proof of Theorem 11.11 

Proof of Theorem \1.1[ Suppose that / € C°° is the solution of 

N 

with compact support. Take the Fourier transform of its both sides. We obtain 

f(o=H(t/\)m/x), (9) 

where = A -1 J2jLo Cje~ 2mdj Z . By Paley and Wiener theorem, /(£) is an 

entire function satisfying /(£) < Cfc|£|~ fc for any positive integer k and £ € R, 
where Ck is a constant. 

By ©, for any M € N, taking the product /(A£)/(A 2 f) • • • J{\ M C)i w e 
deduce that 

M-l 

i/(a m oi = I/col n wt)v (io) 

Suppose that the zero points of H(£) on [0, Dq] are {r\, . . . ,r m }, where Dq is 
the least common multiple of the denominators of dj. Hence, each nonnegative 
root of H(£) has the form of r, + fcDg with some k G Z. By Theorem 13. 1[ we 
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can select a positive number £o such that there exists a c > for which 

UoX 3 - n - kD Q \\ z = U X 3 ~ n\\ z > c 

for every i = 1, . . . , m, and every j = 0, 1, 2, . . .. Therefore, there exists an 
£ > such that H(X j £ Q ) > e a for all j. So, by CGI]), 

|/(A M eo)l > l/te))k M 
Moreover, since / € C°°, for any integer fc, we can find a constant C k such that 

c k (\ M ^r k > i/(A M e )i > im)^ 1 . 

It follows that 

mo)\<C k (X- k ) M ^ k e^ M . 

We can select k so large that A - * < £o- Letting M — > oo, we deduce that 
/(Co) = 0. 

Consider the derivative of /(£) at £o- Note that 
Since /(£o) = 0; we have 

A/'(A£ )=tf(£o)/%). 

It follows that 

M-l 

|A M /'(A A %)| = |/'(£ )| 

Using essentially the same argument, we deduce that /'(Co) = 0. By induction, 
it follows that d J^jP |g=g = for any integer k > 0. Since /(£) is an entire 
function, this yields /(£) = and hence f(x) = 0. This completes the proof. □ 
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4 Proofs of Theorem 12.41 and Proposition 12.51 

To prove Theorem 12.41 we first prove a lemma, which shows that the regularity 
of a refinable function can be determined by a certain smoothness property near 
the endpoints of the support. 

Lemma 4.1. Let f be a compactly supported solution of 

N 

f(x) = Y,c j f(\x-d j ), (11) 
j=o 

with supp f — [A, B], Suppose that there exist a uq £ [A, B] and m € N such 
that f^ m '(x) either does not exist or is discontinuous at uo- Then for any e > 
there are xo £ [A, A + e) and x\ G (B — e, B] such that /' m ^(x) either does not 
exist or is discontinuous at xq and x\. 

Proof. We first consider the case, where xo £ [A, A + e). Without loss of gener- 
ality, we may suppose that do = (see Remark 1). By the result of [B] (which 
was stated in Section 1), we have supp/ c [0, djv(A— l) -1 ]- To prove the lemma, 
note that 

N 

f(X- 1 x)=Y / c J f(x-d J ). 

3=0 

Hence 

N 

f(x) = c - 1 /(A- 1 ^) - cq- 1 Cjf(x - dj). (12) 

Consider the set 

S := {y\ f( m '(x) either does not exist or is discontinuous at y}. 

Since / ^ C m , the set S is nonempty. We choose an element y\ E S. By 
(fl3i|) . {j/i/A,yi — d\, . . . ,2/1 — dpf} n S is nonempty. Therefore, we can find a 
2/2 £ {vi/\ Vi ~ di, ■ ■ ■ , ?/i — cZat} n S such that < j/2 < Vx- By induction, 
there is a sequence yk such that yk € S and < 2/fc < 2/fc-i- If this sequence 
contains 0, then the lemma is proved. Suppose that is not an element of the 



1G 



sequence. By the construction of y k , there exists a fco such that yk a < d\. We 
have 2/fc /A £ S, since yk Q — dj < for 1 < j < N. Thus, for fc > fco, we can 
take j/fc = yk-i/X. It is clear that, for any e > 0, we can find an integer fc so 
large that x :— yk„/X k < £■ 

In case x\ € (B — e, B], we can suppose that d^ — 0. By the same method 
as above, the conclusion follows. □ 

Next, we give the proof of Theorem l2.4l 



Proof of Theorem \2.4\ We suppose that f(x) is smooth on (a,-, ftj+i) where 1 < 
j < M and supp/ C [ai, om+i]- Without loss of generality, we can suppose that 
do = 0, and hence a\ = Q. Let us define 

A k f-(a j )= lim A*/0»0, A k U(aj)= Km A k f(x), 
ax x—Kij- dx ax x->a,j+ dx 

and / fe (aj) = (^) fe /+(a 3 -) - (£) k f-( a j)- We shall prove that either f k { aj ) = 
or lim-c—Kj. {4-) k f{x) exists for all a,j except when fc = fco for some nonnegative 
integer fc . Note that supp/ C [0, om+i]- Then the function /(x) satisfies 

f(x)=cof(Xx-d ) = cof(Xx), xe[0,e] (13) 

for a sufficiently small e > 0. 

We claim that there exists a nonnegative integer fco such that f+°\o) ^ 0. 
For a contradiction, assume that /i_ fc ' (0) = for any nonnegative integer fc. Note 
that /i fe) (0) = for each fc > 0, since f(x) = for x < 0. So / (fc )(aOU=0 exists 
for any nonnegative integer fc. Since fix) is the piecewise smooth function, 
there exists a positive number e\ such that f(x) £ C°°[0, £i). According to 
Lemma |4~TI we have f(x) <E C°°. But then Theorem 11.11 implies that f{x) = 0, 
a contradiction. 

According to (13]), we have /j M (0) = c \ k ° f { + ko) (0). But /j M (0) 7^ 0, so 
c = A~ fc ° . By fig]), we conclude that /(») = \~ k °f(\x) on [0,e]. A simple 
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calculation shows that 

/| fc )(0) = A fc - fc °/j fc) (0) for any k e Z+. 

Clearly, A ^ 1 yields that /| fe) (0) = for k ^ k . Note that /i fe) (0) = for any 
nonnegative integer k. Accordingly, lim x _>o / (%) exists for k ^ fco- By Lemma 
1431 lim x ^ a . /^fa) exists for a11 a 3 if fc 7^ fc o- Put gr(aj) = lim^^. /( feo+1 )(a;) 
and = /( feo+1 ^(a;) for x ^ aj. We will show that g(x) satisfies the following 
refinement equation 

AT 

g(x) = \ k »+ 1 Y,c 1 g(\x-d J ). 

3=0 

Fix .t e R. If {x , Ax - dj,j = 1, . . . , N} n {o Jt j = 1, . . . , Af + 1} = 0, by 
taking the (fco + l)-th derivative at To on both sides of 

N 
3=0 

we obtain 

N 

g(x )=\ ko+1 Y / c 39(^o~d J ). (14) 

3=0 

Let us consider the remaining case when the intersection of two sets is non- 
empty. Without loss of generality, we may suppose that 

{x ,Xxo-dj,j = 1, . . . ,N}D{aj, j = 1, . . . ,M+1} = {x , Xx -dj,j = 1, . . .,JV } 

with an integer N . Select an e > such that 

(x - e, x Q ) n {a jt j = 1, . . . , M + 1} = 

and (At — <4 — Ae, Ato — dfc) n {a-,-, j = 1, . . . , M + 1} = 0, for k = 1, . . . , N. 

Hence f (x) , f (Xx — dj) £ C 00 (xo — e,xo), where j = l,...,N. Then, for 
x € (.to — e, To), we have 

N N 

g{x)-X k ° +l J^c j g{Xx-d j )=X k ° +l J2 WiXx-dj). (15) 

3=0 j=N + l 
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By taking the limits on both sides of (fT5j) . and noting that 

N N 

hm_ X k " +1 c 3 g{Xx - dj) = X ko+1 ^ c jg (\x Q - dj) 

j=N + l j=N + l 

and 

No No 

^imJjW - X k " +1 Cj5(Ax - dj)) = g(x ) - X k « +1 g Cj g(\x - dj), 

3=0 j 

we conclude that 

N 

g(xo)^X ko+1 Y / c 3 9(Xx -d J ). 

3=0 

Combining the results above, we arrive to the equality 

N 

g(x) = A fco+1 ^2 c 39( Xx ~ d j)i f or al1 x e R 

3=0 

Next, by Theorem ll.il we obtain that g(x) = 0, since the function g £ C°° 
is compactly supported. Using f {k " +1) \( a] ,a 3+1 ){x) = g\( aj ,a j+1 )(x) = 0, we 
conclude that f(x) is a spline function. □ 



Now, we begin the proof of Proposition 12.51 For this, we shall give some 
definitions (see also [3J) and a lemma. The functions of the form G(w) — 
^2f=o a,je~ 2 ^ lbjW are referred to as quasi-trigonometric polynomials, where dj £ 
C,dj 5^ and bj £ R such that bo < b\ < ■ ■ ■ < 6 at. In case bj £ Z, such 
polynomials are simply trigonometric polynomials. If &o = 0, G(w) is called a 
normalized quasi-trigonometric polynomial. For the quasi-trigonometric poly- 
nomial G(w), one can write 

G(w) = e~ 2 ' KiriW G l {w) + e~ 2mr2W G 2 (w) + ■■■ + e- 2nir ' w Gi(w), (16) 

where each Gj is a trigonometric polynomial and < r% < • • • < rj < 1 are 
distinct. It is easy to see that up to a permutation of terms this decomposition 
is unique and we shall call (|16[) the standard decomposition of G(w). Moreover, 
as in [3J, we write Ag(w) := X^jLo Cje~ 2mkjW for the greatest common divisor 
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of the trigonometric polynomials {Gj} normalized so that fco = 0, cq = 1 and 
kj > are distinct. It is easy to see that Ac(w) = ae~ 2m: ' w G(w) for some 
constant a and an integer j, if and only if, G(w) is a trigonometric polynomial. 
Then we have 

Lemma 4.2. Suppose that two normalized quasi-trigonometric polynomials G\ 
and G2 satisfy the equality Gi(w)G 2 (w) = R(w)(l — e~ 2nlw ), where R(w) = 
1 + o J jE~ 2 ' nlhjW and < bi < b 2 < ■ ■ ■ < b n Are irrational numbers. Then 

there exists a positive integer P and a set Sq C Si :— {0, 1, . . . , P — 1} such that 
GiW=«(i«)II jeS >" Wi/P -«" 2 ™ /P ) ™dG 2 (w) = p(w)n jeSASo (e-^/ p - 
e~ 27r ™/ p ), where a, /3 are some quasi-trigonometric polynomials. 

Proof. Set G P (w) := R(Pw)(l - e - 2nlPw ), G 1P (w) := G x (Pw) and G 2 p{w) := 
G 2 (Pw) for any integer P. Note that Ag p (w) = 1 — e ~ 27TlPw . We claim that 
the conclusion holds if there exists a positive integer P such that 

A Glp (w)A G2P (w) = 1 - e-w. (17) 

Indeed, on both sides of (fl7|) we have trigonometric polynomials. Setting 
z = e ~ 27rlw , we see that 1 — z p is a product of two polynomials in z. So there 
exists a set S C {0, 1, . . . , P- 1} such that A Glp (w) = UjES (e" 2Try - e^ 2 ™ 1 ), 
which implies G\(w) = II :) es ( e ~ 27r4: '^ P ~ e ~ 27ritu ^ P ) Q! ( u; ): where a(w) is a quasi- 
trigonometric polynomial. Similarly, we have G 2 (w) — rLj<ES 1 \s ( e ~ 27 ™ : '^ P — 
e~ 2niw / p )[3(w), where /3(w) is a quasi-trigonometric polynomial. The claim 
follows immediately. Hence, in order to complete the proof, it suffices to show 
that there exists a positive integer P such that Aq 1p (w)Ag 2P (w) = \ — er 2mPw . 

Write Gi(w) = Y.*l Q a l3 e- 2mb ^ w and G 2 (w) = J^f=o a 2j e- 2 ^ w . Then 
there exists a positive integer, say, P such that each element in the set 

{P6y,P62fc,P(6y - b xj ,),P{b 2k - b w ) | < < N x , < M' < iV 2 } 
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is either integer or irrational. Write 



R(w)(l - e 



-2-iriPw 



) = G XP {w)G 2P {w). 



(18) 



We claim that for each wq satisfying 1 — e 2nlPw <> = 0, one has 



A GlP (wo)A G2P (w ) = 0. 



Suppose that there is an element wo e {w\l — e 2mPw = 0} such that neither 
Aq 1p (wq) nor Ac 2P (wo) is zero. Write wq in the form jo/P + la, where < 
jo < P, j , 1 G Z. Then A Glp (j /P + 1)^0 and A G2P (j /P + 1)^0 for 
every I € Z. 

On the other hand, for each J 6 Z, one has Gip(jo / P + I)G2p(jo / P + 1) = 
(see (|18p). It follows that there exists an infinite set Zo C Z, such that 
Gip(jo/P + I) = for each J e Zo (otherwise, one can replace Gip by G 2 p)- 
According to the choice of P, we have the following decomposition 



where < r± < . . . < r\ are all irrational numbers (except perhaps for r\ = 0) 
and Qj(w) are trigonometric polynomials. Moreover, by the choice of P all the 
differences rj — Tj are also irrational. Substituting jo / P + I with I G Zo into 
Gip(w) we have 



Gip(j /P + I)=J2 e- 2 * irkI e- 2 * irkjo/P Qk(jo/P) = £ e" 2 ^ fc/ T4 = 0, (19) 



where V fc = e~ 2wirkj °/ p Q k (j /P) are all independent oil. Suppose that 7 1( . . . , J; 
are Z distinct elements of Zq. (As Zo is infinite, one can find in it I distinct ele- 
ments.) Let M be the I x I generalized Vandermonde matrix 



Note that all the r^, where k > 2, are irrational and all the differences r k — rj 
are irrational, so the matrix M is non-singular. Taking I — I±, . . . , Ii in (|19|) . 



M = (e 



1-rrirkIj 



l<k,j<h 
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we obtain MV = 0, where V := (Vi, . . . , Vi) T . It follows that all Vj = and 
thus Qk(jo/P) = for every k. Hence Aa lp (j /P + I) = for each I G Z, 
which contradicts to Ag 1p (wo)Ag 2P (wo) ^ and proves the claim. 

Moreover, for each wq such that A<3 lp3 (wo)Ag 2P (wo) — 0, one has 1 — 

e -1^P WQ = q Indccdj jf I _ e -2^Pw„ ^ Q) thcn j _ e -2 7 r J P( Ji . +/) _^ for 

any I £ Z. So, i?(i«o + -0 = for each / G Z. By a similar method, we can 
show that R(w) = 0, which contradicts to the definition of R(w). Similarly, one 
can show that the roots of Ag 1p (w)Ag 2P (w) = are of multiplicity 1. Hence 
we have AG 1P (ti))Ac 2P (ic) = 1 — e - 27 " Ptu . The lemma follows. □ 

Proof of Proposition \2.5\ Suppose B(x\(l, y/h/2)) — fi(x) * f 2 {x), where fi(x) 
and f 2 {x) are splines. By the Fourier transform, we have 

fi(w)f 2 (w) = — . (20) 

Jn > 2iriw 2iriy/5/2w V ' 

By Corollary 2.2 of [3], fj{w) has the form of pj(w)/w, where Pj(w) is a quasi- 
trigonometric polynomial with Pj(0) = for j = 1, 2. Then we have 

Pl (w)p 2 (w)=R(w)(l-e- 2 ™), 

where R(w) := ^—^ V /— — ■ Lemma 14.21 with So (or Si \ So) containing 

— 47r 2 -^/5/2 

implies that there exists an integer Pi such that either pi (w) (or p 2 (w)) is of the 
form of (1 — e~ 2mw ' Pl )qi(w), where qi(w) is a quasi-trigonometrical polynomial. 
Without loss of generality we may assume that pi(w) = (1 — e~ 27Tlw / Pl )qi(w). 
Note that w = is a root of pi(w) and p 2 (w) of multiplicity 1. Hence, by a 
similar method as above applied to w' = y/5/2w, we obtain that P2(w) is of the 
form p 2 (w) = (1 — e ~ 2m VW^ u 'i ' P2 )q2(w) , where P 2 is a positive integer and q 2 
is a quasi-trigonometrical polynomial. 

Set Z(f) := {w\ f(w) = 0, w G C} and Z'(f) := {w\ f(VWw) =0,m£ C}. 
We claim Z(pi) \ Z'(pi) ^ 0, i.e., there exists icqGC such that pi(w ) = but 
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Pi(VlOwo) 7^ 0. To prove this, note that 

Z(p 1 /q 1 )={I 1 P 1 \ h € Z} 

and 

Z'{ Pl q 2 ) = {I2/VT0, I3P2/5 + k/5\ h,heZ,l<k<P 2 - 1}. 

We can see that Z{p l /q 1 ) \ Z'{p 1 q 2 ) ^ 0. Since Z(p 1 /q 1 ) C and Z'{p x ) C 

Z'(pk7 2 ), we derive that Z( Pl ) \ Z'( Pl ) ^ 0, where (^(pi/gi) \ ^'(pi<? 2 )) C 
(Z(p 1 )\Z'(p 1 )). 

Hence there exists a wq such that pi(wo) = while pil^/TOwo) 7^ 0. As a 
result, we see that pi(\/l0w) /p\{w) cannot be a mask polynomial. Thus, f\ is 
not -x/lO-refinable. The proposition follows. □ 



5 Proofs of Theorem 12.61 and Theorem 12.7 



Proof of Theorem \2.b\ Let us begin with part (B). Write A = (mi, . . . ,m„), 
where rrij € R \ 0. We shall prove that, for each element of A, for instance, mi 
there exists an element m £ A and an integer p such that m = pmi/A. 
Consider the Fourier transform 

2 g~ 2-jriwmj 



B(w\A) = ]J 



3=1 3 

Since B{x\A) satisfies the refinement equation, we have 

B(Xw\A) = p(w)B(w\A), (21) 

where p{w) is the mask polynomial. Note that 

B(Xw\A) _ " Q(Am jW ) 



p(w) = £W£ = A -„-Q 



B(w\A) j-Ji Qirrijw) 
where Q(w) = 1 - e r 2mw . Put 

Zj := Qinijiv) =0,ic^0} and := {w\ Q(Xmjw) = 0,w^ 0}. 

23 



Since p(w) is an entire function, one has 

n n 

\Jz 3 c\Jz>. (22) 

A simple calculation shows that 

Z j := {I s / mj \ Ij€Z\ 0} and % := {fcj/(mjA)| % £ Z \ 0}. 

Let us consider Z x . Put J x ■= {j\Z x n 7^ 0}. By we see that Z\ C 

UjGJi ^'r Select an entry in Ji, say, w. Since Z\f\Z' u ^ 0, there exist 2i € Z\0 
and k u G Z \ such that l\jm\ — k u /{m u \). Hence, we can find two coprime 
integers P\ u ,Qiu such that m u = Pi u mi/(Qi u A). Similarly, for any index 
u G Ji, we can find two coprime integers Pi u ,Qiu such that 

™ u = ^ > i«mi/(Qi„A). (23) 

Note that for each fixed I\ G Z \ there exist u G Ji , k u G Z \ such that 

Ii k u k u Qi u 

mi m u X m x P\ u ' 

since Zi C Ujgjj Hence, for any I\ G Z\0, there exist w G Ji and fc„ G Z\0 
such that = P\u/Qiu, i-e., fe u = hPi u /Qi u - So, there exists w G Ji such 

that <5i M |/i for any Ji G Z \ 0. 

We claim that one can find a wi G Ji such that Qi Ul = 1. Indeed, suppose 
Qiu 7^ 1 for each u G Ji- Take Ji = n«eJ + 1- But then /1 is not divisible 

by Qiu for each it G Ji, a contradiction. 

Next, by (|2"3"]) . there exists an index iti such that 

Piu x mi 
m Ut = — - — . 

Setting m := m Ul and p :— P± Ul we complete the proof of (B). 

Consider parts (A) and (C). Using (B), one obtains an infinite sequence 
ui,u 2 ,..., such that m Uk = P UkUkl m Ukl / X, for some P UkUk _ l £ Z, where 
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k > 2. Since each index in this sequence is at most n, the sequence contains 
two equal indices. Without loss of generality, we may suppose that ui+i = u±, 
where / is a positive integer. Then 

^112 V *-U\ U2 m<Ui ) 

^114 ~^Pu3U4™U37 (24) 



This yields A' = P UlU2 • • • -P u ,ui € Z, proving part (A). Moreover, by using 
([24]) , one can see that m„ 2 = m Ul pi/\, m„ 3 = m Ml p 2 /A 2 , . . . , to u , = j3/_im Ml /A i_1 , 
where pi := P Ul u 2 ,P2 := pxP U2 u 3 , ■ ■ ■ ,Pl-i ■= Pl-zPuum-i- Hence |p 3+ i for 
7 = 1,2,...,/ — 2 , and pi_i|A ( . Thus the sub- vector (m Ul , . . . , m^, ) has the form 
of (v, wpiA -1 , . . . , vpi-i A~' +1 ), where w := m Ul , proving (C). 

It remains to prove (D). Consider the Fourier transform of B{x\A) 

n 1 p — 27Tim,jW 

B(w\A) = U—-. . (25) 

- LJ - 2mm jW 

j=i j 

Since B(x\A) is a A-rcfinable spline with integer translations, by Theorem 
12.21 and Theorem 12.31 we can write B(w\A) in the form 

k~ 1 -1 — 2-ni\^w 

B{w\A) = e~ 2 ™° w l[p(e 2 « iX3w )(—— )\ (26) 

j=o 

where zq :— do(l + A + • • • + \ k ~ 1 )/(\ k — 1), h is an integer and p{z) is a 
polynomial. 



Set 



fc-i 

— 2-iri\ J w\h 



3=0 



and 



g 2 (w) = n(i 
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^ — 27TirrijW\ 



Comparing (f2"5)) and (f2"6")l . one gets G2{w) = e 27!law G\{w), where a is a constant. 
Expanding l~[j=o p(e 2nlX3w )(l — e^ 2 ^ 1 ^ w ) h , we see that each term has the form 
^mwT. j=0 bj\ 3 ^. g 2. Since A fc > 1 is an integer, by the theorem of 

Capelli (see [3D] or p. 92 in [M]), A is an algebraic integer of degree k. Thus the 
difference between two distinct numbers in the form of Y2j=o 1S non-integer. 
We conclude that Ag^w) = p(e 2mw )(l — e - 2 ™w^h_ j± s i m ple observation also 
shows that A G2 (w) = T\ mjeZ {l - e - 2mm > w ). Since A Gl (w) = A G2 (w), there 
are h integer entries in the vector A. Without loss of generality, suppose that 
mi, . . . , rrih S Z. Note that 

fc-l fc-l h 

j=0 j=0 r=l 

Then G 2 (w) = e 2mw JJ^Zo lt=i( l ~ e^ 27 "™^')- It follows that a = 0, so 
the vector A can be written as a union of (m r , \m r , . . . , A fe_1 m r ), r = 1, . . . , h, 
proving part (D). □ 

Proof of Theorem \2.7\ We begin with (A). Assume that 

M = (mi,...,m„) G R sxn . 
It is well-known that the Fourier transform of B{x\M) is 

^ 2iriw T mj 



B{w\M) = Yl 



2iriw T mj 
j=i J 

Since B{x\M) is A-refinable, the mask polynomial H is given by the formulae 
H(w) = B(Xw\M)/B(w\M). 

Put W :— {w\w T rrij ^ 0, for all j}. Select a, wq G W and consider 

n i g — 27rizw^ rrij 

f(z\w ) := B(zw \M) = [T , T , (27) 

where zeR. Observe that /(z|wo) can be considered as the Fourier transform 
of the univariate box spline B{x\{w^mi : . . . ,ui^m n )), which is A-refinable for 
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each fixed wo £ W with the function f(\z\wo)/ f(z\wo) = H(zwq) being a quasi- 
trigonometric polynomial. Then, by Theorem l2.61 there exists a positive integer 
I such that A' £ Z, proving part (A). 

Select an entry in the vector (wqUIi, . . . , w^m„), for instance, w^rriQ. Since 
B{x\(wq m\, . . . , WQin n )) is A-rehnable, one can find an integer p Wo and an index 
j Wo satisfying 2 < j Wo < n such that w r mj mg = p Wo w^mo/\ for any wq £ W 
(see part (B) in Theorem 12. 6[) . Each wo corresponds to an index j Wo . For an 
index j, set W(j) = {w Q £ W\j Wo = j}. Then Uj =2 W(j) = W. Since the s 
dimensional Lebesgue outer measure of W is infinite, there exists a subset of 
W, say, Wq such that, for any wq £ Wq, the index j Wo is a constant and the s 
dimensional Lebesgue outer measure of Wo is positive. We suppose the constant 
index j Wo is u\. Then 

WqITIo rp 

Pw — t — = w m ui (28) 
for any wo £ Wo. We now consider all integers p Wo . For each q £ Z, set 

W Q (q) := {w \p Wa = q, wo £ Wq}. 

Then 

(J W (q) = W . 

qez 

We claim that there exists a positive integer q\ such that the s dimensional 
Lebesgue outer measure of Wo{qi) is positive. Indeed, the s dimensional Lebesgue 
outer measure of Wq is positive and the set Z is countable. Hence there exist 
s linearly independent elements of Wo(qi). Let us denote them by wi, . . . ,w s . 
We have 

wjmo _ 
9i — ^ — = Wj m Ul , 

for each 1 < j < s, by ([28]) . It follows that 



wJV = 0, (29) 
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for each 1 < j < s where V :— (qi^ — m Ul ). Let A be the s x s matrix 

A = [wj)x<j< a . 

Then A is non-singular, because Wj are linearly independent. The equality (|29p 
can be written as AV = 0. It follows that V = 0, because A is non-singular. 
Hence m Ul = qimo/X. Putting m := m Ul ,p :— qi, we complete the proof of 
(B). 

Part (C) can be proved by the same method as in the proof of the part (C) 
of Theorem 12.61 We omit the details. 

It remains to prove (D). Set W := W H Z s . By the definition of W, one 
can see that W = Z s \ (U^ify), where H s := {w £ Z s \w T mj = 0}. Then 
B(x\(wQirii, . . . ,WQm n )) is A-refinable with integer translations for every wq £ 
W. By Theorem 12. 61 the vector (wqITIi, . . . ,WQm n ) can be written as 

rp T T T T T \ 

[Wq m Ull , . . . , Wq rn Ulk , w m U21 , . . . , w rn U2k , . . . , w m Utl , . . . , w m Utk ). 

Here, for each fixed r = 1, . . . ,i, one has w r m Urh+1 = Awq m Ur k , where h = 
1, . . . , k — 1 and t :— n/k is an integer. Hence, each wq £ W corresponds to the 
index vector P(wq) := (itn, . . . , Uifc, . . . , . . . , u t k), which is a permutation 
of (1, ...,n). 

Let us denote the set consisting of all permutations of (1, . . . , n) by P. For 
each p £ P put 

W'(p) :={w Q £ W'\P(w )=p}. 

We claim that there exists a po £ P such that span(W / '(po)) = R s , i-e., that 
there are s linearly independent vectors in W (po)- For a contradiction, assume 
that span(W(po)) is contained in a (s — l)-hyperplane for every po £ P. Note 
that 

U^pW'ip) = W' = Z S \ (U] =1 H 3 ). 



2N 



Then 



Z s = (U peP f'(p)) U (Uj^Hj). 



(30) 



Since #P is finite, the equation ([50)) shows that Z s can be written as a finite 
union of hyperplanes, yielding a contradiction. 

Without loss of generality, we may suppose that span(W'(po)) = R s for 
Po = (tin, . • • , wife, . . . , Uti, . • • , ittfe)- Then, one can select s linearly independent 
vectors in W(po), say wi,...,w s , such that for each fixed 1 < r < t and 
1 < j < s, one has wJm Ur h+1 = XwJm Urh , where 1 < /i < fe — 1. Hence, for 
fixed r and h, we obtain the following linear equations 



where j — 1, . . . , s. Solving these linear equations, we get m Ur h+1 = Xm Ur h , 
where r = 1 , . . . , t and h = 1, . . . ,k — 1 . 

Hence the matrices (m Url , . . . , m Urk ), 1 < r < t, are of the form (mo, Amo, . . . , 
A fe_1 mo) with mg := m Url . Consequently, the matrix M can be written as a 
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